Chiral condensate with topological degeneracy in graphene and its manifestation in edge states 
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Role of chiral symmetry in many-body states of graphene in strong magnetic fields is theoretically studied 
with the honeycomb lattice model. For a spin-split Landau level where the leading electron-electron interaction 
is the nearest-neighbor repulsion, a chiral condensate is shown to be, within the subspace of ;i = Landau level, 
an exact many-body ground state with a finite gap, for which calculation of Chem numbers reveals that the 
ground state is a Hall insulator with a topological degeneracy of two. The topological nature of the ground state 
is shown to manifest itself as a Kekulean bond order along armchair edges, while the pattern melts in the bulk 
due to quantum fluctuations. The whole story can be regarded as a realization of the bulk-edge correspondence 
peculiar to the chiral symmetry. 
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Introduction. — While the physics of graphene started from 
the one-body electronic structure as a Dirac fermion, a pos- 
sible relevance of electron correlation in graphene has been 
intensively studied after a gap opening in the n - Q Lan- 
dau level (LL) was experimentally observed in strong mag- 
netic fieldsii2i Since it is difficult to explain the gap within 
a simple one-body problem, considerable theoretical efforts 
have ensued to clarify many-body effects in graphene quan- 
tum Hall (QH) systems.—— However, little attention has been 
paid on how many-body effects should reflect the chiral sym- 
metry in the graphene QH regime, which is after all a fun- 
damental symmetry inherent in graphene's honeycomb lat- 
tice. On the one-body level, the effects of chiral symmetry 
in graphene is well understood: To start with, the symme- 
try guarantees the emergence of doubled Dirac cones in the 
Brillouin zone, which can be interpreted as a two-dimensional 
analogue of the Nielsen-Ninomiya theorem well-known in the 
four-dimensional lattice gauge theory. We can even exam- 
ine the wave functions in terms of Aharonov-Casher argu- 
ment, which states that chiral symmetry topologically pro- 
tects the degeneracy of the « = LL against random gauge 
fieldsJS A similar situation occurs for ripples in a graphene 
sheet, which can be modelled by random hopping amplitudes. 
Kawarabayashi et al. have shown that the n - Q LL exhibits 
an anomalously sharp (delta-function-like) density of states 
(DOS) as soon as the spatial wavelength of the ripple exceed 
a few lattice constantsJi 

In the presence of electron-electron interactions, on the 
other hand, the role of chiral symmetry has primarily been 
investigated in zero magnetic fields in the context of sponta- 
neous symmetry breakingAi^r— While these studies mainly 
employ a Dirac field model in a continuum space to dis- 
cuss many-body gap formation, such an effective treatment 
may well overlook the essence of graphene's chiral symme- 
try, which is intimately related to the underlying honeycomb 
lattice. 

With this background, we shed light in the present work on 
how the chiral symmetry influences the many-body problem 
in graphene QH effect, by fully taking account of the lattice 
structure. We first examine the many-body problem with ex- 



act diagonalization in a subspace projected onto the n - Q LL. 
Working on the subspace enables us to classify many-body 
states according to a notion of the total chirality of the filled 
zero modes. In terms of this, for a "bipartite" electron-electron 
interaction such as the nearest neighbor repulsion, which is 
the dominant interaction for a spin-split LL, we show that the 
many -body ground state is exactly identified to be a chiral con- 
densate with a topological degeneracy of two. We confirm 
numerically that there exists a finite energy gap to the first- 
excited state, which makes the Chem number of the ground 
state well-defined. The total Chern number contributed by the 
filled zero modes along with the negative energy states ("Dirac 
sea") turns out to be zero, which implies the system is a Hall 
insulator with vanishing Hall conductance. 

Despite the cancellation of the Chern number in the bulk, 
however, we move on to show that the topological nature of 
the chiral condensate is in fact made manifest as an emer- 
gence of a Kekulean bond order in the edge state along arm- 
chair edges of the honeycomb lattice (in sharp contrast to 
zigzag edge states in one-body problem). On the mean-field 
(MF) level the Kekule pattern is shown to appear in the bulk 
with a Kekulean degeneracy of three, so that the present re- 
sult amounts that the mean-field order, while dissolved in the 
topologically-degenerate chiral condensate in the bulk, resur- 
faces along an armchair edge. This can be interpreted as an 
example of bulk-edge correspondence,- which states that a 
topologically nontrivial bulk state should always accompany 
a characteristic edge state. In the case of graphene, the bulk- 
edge correspondence becomes peculiar in that the chiral con- 
densate manifests itself as a freezing of a Kekulean bond order 
along a specific (armchair) edges. 

Chiral symmetry. — To model interacting electrons on a 
honeycomb lattice, we assume that spin degeneracy is Ufted 
by a Zeeman splitting, so that the leading Coulomb interac- 
tion reduces to the nearest-neighbor repulsion. The Hamil- 
tonian then reads l-i - 'H^n + 'J^im, where the kinetic term 
'^kin = -tYjiifiis'^^'-'clcj + H.c.) s c'"//kinC dcscribcs electron 
hopping between adjacent sites (ij) with strength f > 0. The 
magnetic field is included as the Peierls phase Oij such that 
magnetic flux per elementary hexagon equals 2q 6*,^ - litip 
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in units of a magnetic flux quantum hje. For a honeycomb 
lattice with A^.(o) sites in sublattice •(o), c' = (c',,c[) with 
c^(^j a row of creation operators for sublattice "(o) and //kin 
is a square matrix of dimension A^, + A^o- If we introduce 
r = diag(/,,-/o) with an identity matrix /,(o) of dimension 
A^.(o), the kinetic term satisfies an anticommutation relation 
{//kin, r) = 0, which defines the chiral symmetry. The symme- 
try implies that, if is the A:-th eigenvector with energy e^, a 
chiral partner Yxj/k exists with an energy -e^ . This makes the 
n = LL special in that we can take (/^^ as an eigenstate of F 
as T\l/k± - +i//k±- The interaction between spin-polarized elec- 
trons is expressed in a particle-hole symmetric form as 'Hint - 

const., where Vij is the strength of electron-electron interac- 
tion, n, s c'. Cj the number operator at site /. 

Chiral condensate. — We start with an investigation of the 
many-body problem at half filling. Since it is difficult to treat 
all the many-body states exactly, we shrink the Hilbert space 
by projecting onto the « = LL. Such a treatment is valid as 
long as \ Vij\ is perturbatively small compared with the Landau 
gaps around the n - Q LL4^ In the « = LL, we take a zero 
mode multiplet if/ = where we have decomposed it 

into eigenstates of the chiral operator, iff± = (iAi±> ' ' ' ^4'Mi±) 
with degeneracy M±. Note that the zero modes are localized 
on each of the sublattices as = -^^ ( ) and = ^ ° 
If we introduce the negative-energy multiplet (p - {(fi,ip2,- ■ ■) 
such that Hi^mifk - ekfk with < 0, (if/+,if/-,ip,rip) form 
a complete set, so that the fermion operator is expanded as 
c - i/zd + ipd< + Tipd>, where = ( ) and d^ are columns of 
fermion operators for the zero modes and the e S states, re- 
spectively. In the projected subspace, the total Hamiltonian 
is written, up to a constant, as ■// = \Tjiti^iA^'i^'j^fj + 
CiCjC^.c\), where the projected fermion operator, c, = ij/d, no 
longer satisfies the canonical anticommutation relations. Still, 
the generator of the total chirality of the filled zero modes can 
be defined as ^ = c'Tc - d\d+ - did-. Due to the invariance 
of 'H for the chiral transformation H 1-^ e'^^iie-'"^ = "H, 
the total chirality is conserved, YH, 0] - 0. Then the many- 
body states are classified according to the total chirality Xtot- 
For repulsive interactions Vjj > 0, 'K is semi-positive def- 
inite. Furthermore, when the interactions are bipartite (i.e., 
only act between different sublattices, Vie,.je» = Vieojeo - 0) 
as the case with the nearest-neighbor interaction, chiral con- 
densates |G±) = d'l^ ■ ■■d'ljD,) (with |D<) s n„,(cV),«IO) 
denoting the Dirac sea) constitute a ground state doublet 
T = (|G+), |G-)), since c,c,|G±) = c';c"i"|G±) = for / e • and 

j 6 o.i^il We call this a topological degeneracy of two. Un- 
like a simple charge density wave (CDW), one may mix |G+) 
and |G_) through a unitary transformation i-> ^"^(d with 
cl) e U(2) even in a finite system. One can numerically con- 
firm that the chiral condensate remains the ground state unless 
the non-bipartite potential is large,— even though the ground- 
state energy is nonzero. In the rest of the paper, we focus on 
the leading interaction, i.e., the nearest-neighbor repulsion for 
simplicity. 

Many-body gap. — We next calculate excitation energies nu- 




FIG. 1 : (color online) Magnetic-field dependence of the many-body 
gap A in the n = LL. The result is displayed for 30 electrons with 
the consecutive values of the linear dimension of the system, L, con- 
nected by a line with the results for L = 3 (mod 3) forming a lower 
envelope (A oc 0^). The inset shows the gap at L = 3/ against 1/L 
for the magnetic flux (p = 1/12 (circle), = 1/27 (triangle), and 
= 1/48 (square). 



merically with the exact diagonalization method. In the pro- 
jected subspace, the strength of the nearest-neighbor repulsion 
y > is the only energy scale, which acts as the unit of en- 
ergy. Full energy spectra of '7Y for finite systems suggest that 
the first excited state appears in the sector of ;^ftot = ±{M+ - 2), 
which is created from a chiral condensate |G±) by single 
chirality-flippings analogous to the projected single-mode ap- 
proximation. "^""^ Noticing this, we further shrink the Hilbert 
space by focusing on the sector of Xtot = ±{M± - 2). This 
enables us to obtain the energy of the first excited state, or the 
energy gap A, with calculation cost of the order of 0{M^). We 
consider a system on a torus composed of 2L^ lattice sites with 
a linear dimension L. For investigating a weak-field regime, 
we adopt the string gauge^2i where a magnetic flux is given 
by <p = m/L- with an integer m (> 0) and M+ - m zero modes 
are obtained for each chirality. 

In Fig.[Tl we plot A for 30 electrons as a function of 4> with 
L changed consecutively. We immediately notice that the re- 
sult exhibits a marked periodicity of three, where the values 
for L = 2i (mod 3) form a clear lower envelope with a scaling 
A oc (p^, while those for L 31 deviates from this. The latter 
behavior is considered to be a finite-size effect, since the devi- 
ation diminishes with the sample size. To confirm the scaling, 
the inset plots A/0^ at L = 3/ against 1 /L for various values of 
(/>, which indicates the scaling law A oc (/i^ is very accurately 
obeyed. 

Hall conductance. — Let us now consider Hall conductance 
of the chiral condensate. By the Niu-Thouless-Wu formula.jS 
the Hall conductance is written with the Chern number— as 

h No 2ni J 

where No is the degeneracy and A is the non-Abelian Berry 
connection that describes multiplets.-'* In terms of the ba- 
sis that diagonalizes 0, we have C = C+ + C- with 
C± - J{dG±\dG±). Each term is further decomposed 
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FIG. 2: (color online) Mean-field results for the energy spectrum 
(a), its blowup around the n = LL (b), and the bond strength |A,j| 
plotted in a real-space (c). The parameters are L = 15, (p = 1/15, and 
V/t = 0.25. A dashed hexagon in (c) is an enlarged unit cell with 
arrows primitive vectors. 

as C± - C^^ + Cd^ with C^^ - J Tr M^dt//ldt//± and 
Cd^ - 2^ /Trii(/j' ii(^. By the charge conjugation, we have 
+ cZ = -(Q- + CdJ with Cd, = /Tr(r#)'Td^ = 
Co^ ■ Thus the total Chern number of the ground-state dou- 
blet vanishes as C = C^^,^ + + 2Co^ = 0, which may be 
called a topological cancellation. This implies that the chi- 
ral condensate is a Hall insulator with a nontrivial topological 
degeneracy No - 2. 

Bond order. — As have been confirmed in various systems, 
while topological phases are featureless in a bulk, they show 
characteristic boundary states^i^ So a natural question we can 
pose here is: do the edge states in the present system exhibit 
special features despite the bulk Chem number being zero? 
Before presenting the result, however, let us first have a look 
at the mean-field state in the present system in the bulk, which 
will turn out to be instructive. One virtue of a mean-field pic- 
ture2^ is that we can introduce a bond order. Ay = y(cj cj), for 
adjacent sites The dominant part of the MF Hamilto- 

nianis given by TYmf = - Zi</7>[(fe'''" -i-A*pCj!cy-i-H.c.], where 
Aij is determined self-consistently by diagonalizing 'Kmf- A 
spontaneous symmetry breaking is induced by the many -body 
effect for weak magnetic fields, where the density of states has 
a sharp peak at the Fermi energy. In Fig. |2] we show a typi- 
cal MF result for the ordered phase. The energy spectrum is 
plotted in panel (a), where the qualitative structure of the LLs 
is preserved. This comes from the fact that the convergent or- 
der parameters turn out to retain the initial Peierls phase as 
Aij - \Aij\e^'^'-i . The influence of the electron-electron interac- 
tion appears most prominently in the n - LL, where a finite 
gap of the order of opens as shown in the blowup Fig.|2lb). 
To see how the symmetry is broken in the mean field, we show 
in panel (c) a real-space image of the bond order | A,y|, which is 
seen to exhibit a Kekule pattem.i^Si^ This makes the unit cell 
enlarged, which causes K and K' points to be coupled, and 
this in turn opens a finite gap. In this sense we can regard this 
a Peierls transition in the honeycomb lattice. 

On the other hand, the chiral condensate with its topological 
degeneracy of two exhibits in the bulk no bond order as we 



(a) armchair edge |A,y| / V (b) zigzag edge |Ay| / V 




FIG. 3: (color online) Bond strength near an armchair (left) or zigzag 
(right) edge of the doubly-degenerate chiral condensate. Magnetic 
flux is if) = 1/192 for which the magnetic length is /g ^ 8.9a. The 
bond order decays in a bulk away from the edges. 



have seen in Fig. |2] Due to the quantum fluctuation, the bond 
order of the mean field is destroyed and the quantum liquid 
ground state is realized. 

Edge and defect states. — We are now in a position to ask 
the question: what kind of edge states does the chiral conden- 
sate accommodate? Based on the bulk-edge coiTespondence, 
we may expect a non-trivial behavior of the many-body states 
near edges. A prime example is the fractional QH states in a 
2DEG, where a CDW-like behavior emerges along edges of 
a ribbon,^" while in the bulk it melts into the Laughlin liquid 
with no long-range order, which has a q-fo\d degeneracy of 
the fractional QH states at filling v - l/q. Note that a honey- 
comb lattice with edges has QH edge states whose mode lies 
in a LL gap. To perform the projection onto the n - LL 
we set an energy cutoff, the choice of which is shown to have 
little influence on the edge states shown below. 

In Fig. [3] we show |A,y| for the chiral condensate plotted in 
a real space near armchair and zigzag edges. In panel (a), 
we can see that a Kekule-type bond order reminiscent of the 
mean-field result in Fig. |2] emerges along the armchair edge. 
This is the key result in the present work. The enhancement 
in bonds rapidly decays away from the edge in a few lattice 
constants, and |A,y| slightly oscillates with a length scale of 
the order of the magnetic length Ib ~ a/ with a being the 
interatomic spacing. This may naively seem to be analogous 
to the fractional QH edge states in a 2DEG, but here the hon- 
eycomb lattice structure is essential in the ground state. In- 
deed, the ring pattern is locked along the armchair edge in a 
Kekule pattern, while this is not the case with zigzag edges 
[see panel (b)]. In the latter case, the ring pattern is blurred 
by the translational symmetry along a zigzag edge, and a very 
weak stripe pattern parallel to the edge appears. These pat- 
terns related with the three-fold degeneracy of the Kekule pat- 
tern are washed out in the bulk chiral condensate. All these 
are a specific property of a honeycomb lattice model. 

We can further endorse that the lattice structure is at the 
core by looking at the states around lattice defects. When a 
single atom is removed from the bulk honeycomb lattice, one- 
body localized zero modes appear that are protected by the 
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FIG. 4: (color online) Bond strength of the chiral condensate near a 
divacancy composed of two adjacent sites missing. Magnetic flux is 
(p = 1/1200 for which the magnetic length is Ig - 22.3a. The bond 
order reflects twofold axial symmetry of the divacancy. 

chiral symmetry ^^i^ In the presence of the electron-electron 
interactions, however, local chiral symmetry breaking occurs 
spontaneously to lower the energy by inducing effective hop- 
ping in the same sublatticei^i^ Then what if two point defects 
come close to each other? Such a divacancy consists of two 
adjacent missing atoms, and is recently observed experimen- 
tally in ion-irradiated carbon samplesj^ We expect that the 
chiral symmetry may be partially recovered with a reconfig- 
uring of the two symmetry-breaking bonds. We plot in Fig. |4] 
the bond order for the chiral condensate near the divacancy. 
We do confirm that enhancement of the bond order near the 
divacancy which can be considered due to the revival of the 



chiral symmetry. 

We have thus shown that the bond order emerges along 
edges and around vacancies, despite the topological cancel- 
lation C = that might first seem to wipe out any signature 
of the chiral condensate. Indeed, the charge density {«,) itself 
is uniform for the chiral condensate even along edges, which 
is due to the invariance of the chiral condensate for the charge 
conjugation. Thus it is the bond order |A,j| that we have to 
look at as a probe for the chiral condensate. Thus the bond 
order provides a new probe for the many-body effect in half- 
filled graphene applied a magnetic field. The bond order near 
edges should be observable experimentally with some imag- 
ing techniques such as Green's function scanning tunneling 
microscopes?22i^ Since the amplitude of |A,j|/y = Kcjcj)] is 
of the order of the magnetic flux 0, the magnetic field should 
have significant magnitudes. 

Summary. — The many-body ground state at half filling in 
the honeycomb lattice is identified as a doubly-degenerate chi- 
ral condensate for a spin-split Landau level. The many-body 
effect opens a finite energy gap, which makes the chiral con- 
densate a generic topological insulator However, the system 
has a peculiar manifestation of the bulk-edge correspondence 
in topological systems as an emergence of a bond order with a 
Kekule pattern along armchair edges in an exact ground state, 
while the pattern is dissolved in the bulk. 
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